Symmetry, Integrability and Geometry: Methods and Applications 



SIGMA 3 (2007), 060, 31 pages 



Generating Operator of XXX or Gaudin Transfer 
Matrices Has Quasi-Exponential KernelE 

Evgeny MUKHIN 't, Vitaly TARASOV n and Alexander VARCHENKO^ 

t Department of Mathematical Sciences, Indiana University - Purdue University Indianapolis, 
402 North Blackford St, Indianapolis, IN 46202-3216, USA 



E-mail: mukhin@math.iupui.edu, vtarasov@math.iupui.edu 

o : 

•f St. Petersburg Branch of Steklov Mathematical Institute, 



Fontanka 27, St. Petersburg, 191023, Russia 



E-mail: vt@pd mi.ras.ru 



| & Department of Mathematics, University of North Carolina at Chapel Hill, 

ITi ; Chapel Hill, NC 27599-3250, USA 

E-mail: \anv@email.unc.edv\ 

Received March 28, 2007; Published online April 25, 2007 



Original article is available at http://www.emis.de/journals/SIGMA/2007/060/ 



< 

| Abstract. Let M be the tensor product of finite-dimensional polynomial evaluation Y(qI n )- 

■ N 

modules. Consider the universal difference operator 35 = Yl (—l) k 1k{u)e~ kdu whose coef- 

fc=0 

ficients 1k( u ) '■ M — > M are the XXX transfer matrices associated with M. We show 
that the difference equation Df = for an ill- valued function / has a basis of solutions 
consisting of quasi-exponentials. We prove the same for the universal differential operator 



N 

D = ( — l) k Sk(u)du~ k whose coefficients Sk(u) : M — > M. are the Gaudin transfer 

fe=0 

matrices associated with the tensor product M. of finite-dimensional polynomial evaluation 
$l N [xj-modulcs. 
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1 Introduction 



Tb i/ie memory of Vadim Kuznetsov 



In a quantum integrable model one constructs a collection of one-parameter families of commu- 
ting linear operators (called transfer matrices) acting on a finite-dimensional vector space. In this 
paper we consider the vector-valued differential or difference linear operator whose coefficients 
are the transfer matrices. This operator is called the universal differential or difference operator. 
For the XXX and Gaudin type models we show that the kernel of the universal operator is 
generated by quasi-exponentials or sometimes just polynomials. This statement establishes 
a relationship between these quantum integrable models and that part of algebraic geometry, 
which studies the finite-dimensional spaces of quasi-exponentials or polynomials, in particular 
with Schubert calculus. 

We plan to develop this relationship in subsequent papers. An example of an application of 
this relationship see in [5]. 



*This paper is a contribution to the Vadim Kuznetsov Memorial Issue 'Integrable Systems and Related Topics'. 
The full collection is available at |http://www.emis.de/journals/SIGMA/kuznetsov.html| 
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Consider the complex Lie algebra qI n with standard generators e ab , a, b = 1, . . . , N. Let A = 
(A^ 1 ), . . . , A^ n )) be a collection of integral dominant gljy-weig hts, where AW = (Ap , . . . , A$) 

for i = 1, . . . , n and A^ € Z^o- Let M = M A (i) (g> ■ ■ ■ (gi M A ( n ) be the tensor product of the 
corresponding finite-dimensional highest weight gl^y-modules. Let K = (K ab ) be an JV x JV 
matrix with complex entries. Let z\, . . . , z n be distinct complex numbers. For a,b = 1, . . . , N, 
set 

X ab (u, d u ) = 5 ab d u -K ab -) ba , 
where <9 M = d/du. Following |13| . introduce the differential operator 

T>K,M( u idu) = ^ ( _1 ) <T X lai (u,d u ) X 2 a 2 (u,d u ) ••• XNcr N (u,du). 

o-eSjv 

The operator VK,M(u,d u ) acts on M-valued functions in u. The operator is called the universal 
differential operator associated with M, K and z\, . . . ,z n . 
Introduce the coefficients Sq,k,m{ u )-> ■ ■ ■ ,<Sn,k,m(u): 

N 
k=0 

The coefficients are called the transfer matrices of the Gaudin model associated with M, K and 
z\, . . . , z n . The transfer matrices form a commutative family: [S$;^k,m{ u )i Si^k,m{v)} = for all 
k, I, u, v. 

A quasi-exponential is a finite sum of functions of the form e Xu p(u), where p(u) is a polyno- 
mial. 

Theorem 1.1. The kernel of the universal differential operator T>k,m(u,9 u ) is generated by 
M-valued quasi- exponentials. Possible exponents A appearing in the kernel are eigenvalues of 
the matrix K. In particular, if K = 0, then the kernel is generated by M-valued polynomials. 

For K = 0, the theorem was conjectured in [1] and proved in [5]. 

We describe possible degrees of polynomials p(u) appearing in the kernel further in the paper. 

The universal differential operator has singular points at z\, . . . , z n . We describe behavior of 
elements of the kernel at these points. 

The tensor product M may be naturally regarded as the tensor product of polynomial evalua- 
tion modules over the current algebra flljvM- Then the operators X ab (u), T>K,M(u,d u ) may be 
naturally defined in terms of the 0[jv[x]-action on M. 

Similarly, the tensor product M may be naturally regarded as the tensor product of polyno- 
mial evaluation modules over the Yangian Y(gl N ). Then one may define a linear iV-th order 
difference operator acting on M-valued functions in u. That operator is called the universal 
difference operator. Its coefficients commute and are called the transfer matrices of the associa- 
ted XXX model. 

We prove that the kernel of the universal difference operator is generated by quasi-expo- 
nentials. We describe the quasi-exponentials entering the kernel and their behavior at singular 
points of the universal difference operator. 

The paper has the following structure. In Section [2] we make general remarks on quasi- 
exponentials. 

In Section [3] we collect basic facts about the Yangian Y(qI n ), the fundamental difference 
operator and the XXX transfer matrices. In Section 13.41 we formulate Theorem 13.81 which 
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states that the kernel of the universal difference operator is generated by quasi-exponentials. 
Theorem 13.81 is our first main result for the XXX type models. 

In Section [J] we prove a continuity principle for difference operators with quasi-exponen- 
tial kernel. Under certain conditions we show that if a family of difference operators has 
a limiting difference operator, and if the kernel of each operator in the family is generated 
by quasi-exponentials, then the kernel of the limiting difference operator is generated by quasi- 
exponentials too. 

Section [5] is devoted to the Bethe ansatz method for the XXX type models. Using the Bethe 
ansatz method we prove the special case of Theorem 13.81 in which M is the tensor product of 
vector representations of qI n . Then the functoriality properties of the fundamental difference 
operator and our continuity principle allow us to deduce the general case of Theorem 13.81 from 
the special one. 

In Section [6] we give a formula comparing the kernels of universal difference operators asso- 
ciated respectively with the tensor products M AW <g> ■ ■ ■ (g> M A(n) and M A (iw ai ) ® • • • <g) M A („)( a ^, 

where oi, . . . , a n are non- negative integers and AW (aj) = (A^ + Oj, . . . , Ajy + a^) . 

In Section [7] we describe the quasi-exponentials entering the kernel of the universal difference 
operator and the behavior of functions of the kernel at singular points of the universal difference 
operator. Theorems 17.11 17.21 and 17.31 form our second main result for the XXX type models. 

In Sections [HHT21 we develop an analogous theory for the universal differential operator of the 
Gaudin type models. 

Section [8] contains basic facts about the current algebra gt^x], the universal differential 
operator and the Gaudin transfer matrices. We formulate Theorem I8.4| which states that the 
kernel of the universal differential operator is generated by quasi-exponentials. This theorem is 
our first main result for the Gaudin type models. 

Section[9]contains the continuity principle for the differential operators with quasi-exponential 
kernel. 

Section [TU] is devoted to the Bethe ansatz method for the Gaudin type model. 

In Section[TT]we compare kernels of the universal differential operators associated respectively 
with tensor products M A(1) ® • • • fg) M A{n) and M AW ^ <g> • • • ® M A ( n ), s. 

In Section [12] we describe the quasi-exponentials entering the kernel of the universal differ- 
ential operator and the behavior of functions of the kernel at singular points of the universal 
differential operator. Theorems 112.11 T12.2I and 112.31 form our second main result for the Gaudin 
type models. 

2 Spaces of quasi-exponentials 
2.1 Quasi-exponentials 

2.1.1. Define the operator r acting on functions of u as (r/)(n) = f(u + 1). A function f(u) 
will be called one-periodic if f(u + 1) = f(u). Meromorphic one-periodic functions form a field 
with respect to addition and multiplication. 

2.1.2. Let Q be a nonzero complex number with fixed argument. Set Q u = e uln ® . We have 
tQ u = Q U Q. 

Let p € C[u] be a polynomial. The function Q u p will be called a (scalar) elementary quasi- 
exponential in u. A finite sum of elementary quasi-exponentials will be called a (scalar) quasi- 
exponential. 

Let V be a complex vector space of finite dimension d. A V -valued quasi- exponential is a 
T^-valued function of the form ^ a f a (u)v a , where f a (u) are scalar quasi-exponentials, v a G V, 
and the sum is finite. 
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We say that a quasi-exponential Yl,ab Qa u v ab is of degree less than k if v a b = for all b ^ k. 

2.1.3. For given End (V)-valued rational functions Aq («),..., A^{u) consider the difference 
operator 

N 

V = ^A k (u)r~ k (2.1) 

fc=0 

acting on V- valued functions in u. 

We say that the kernel of D is generated by quasi- exponentials if there exist Nd quasi- 
exponential functions with values in V such that each of these function belongs to the kernel 
of D and these functions generate an Ati-dimensional vector space over the field of one-periodic 
meromorphic functions. 

The following simple observation is useful. 

Lemma 2.1. Assume that a quasi- exponential ^2 ab QaU b v ab , with all numbers Q a being different, 
lies in the kernel of® defined in (|2.ip . Then for every a, the quasi- exponential u b v ab lies 

in the kernel ofD. 

The lemma follows from the fact that exponential functions with different exponents are 
linearly independent over the field of rational functions in u. 

3 Generating operator of the XXX transfer matrices 
3.1 Yangian Y(gl N ) 

3.1.1. Let e ab , a,b = 1,...,N, be the standard generators of the complex Lie algebra gl N . We 
have gl N = n + © () © n~, where 

n + = ffi a <fe C • e ab , f) = ®% =1 C • e aa , rT = ® a>b C • e ab . 

For an integral dominant q (^-weight A € fj*, denote by M\ the irreducible finite dimensional 
gljy-module with highest weight A. 

For a gl^-module M and a weight fi 6 fj*, denote by M \p] C M the vector subspace of 
vectors of weight [i. 

3.1.2. The Yangian Y(gl N ) is the unital associative algebra with generators , a, b = 1, . . . , N 
and s = 1, 2, . . . . Let 

oo 

T ab = Sab + T iS~ S > a,b=l,...,N. 

s=l 

The defining relations in Y(qI n ) have the form 

(u - v) [T ab (u),T cd (v)\ = T cb (v)T ad {u) - T cb (u)T ad (v) , 
for all a, b, c, d. The Yangian is a Hopf algebra with coproduct 

N 

A : T ab (u) h-> y^T cb (u) <S> T ac {u) 

c=l 

for all a, b. 

3.1.3. We identify the elements of End (C^) with N x iV-matrices. Let E ab G End (C^) denote 
the matrix with the only nonzero entry 1 at the intersection of the a-th row and 6-th column. 
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Let P = Ea,6 E ab 8> E ba , R(u) = u + P £ End (C N ® C*) and T(u) = £ o6 E ab ® T ab (u) € 
End(C Ar ) ® y(g[^y)((u _1 )). Then the defining relations for the Yangian can be written as the 
following equation of series in u - 1 with coefficients in End (C^) ® End (C^) ®Y(gt N ), 

R (m {u _ y) T (13) (n) r (23) {y) = r (23) (u) T (13) {u) R (12) ( 

3.1.4. A series f(u) in u -1 is called monic if f(u) = 1 + 0(u~~ l ). 
For a monic series /(it), there is an automorphism 

Xf : Y( Q l N ) - y(gl^), T(u) ^ /(«)T(u). 

The fixed point subalgebra in y(jj[jy) with respect to all automorphisms Xf is called the Yangian 
Y(sItv). Denote by ZY(gl N ) the center of the Yangian Y(gl N ). 

Proposition 3.1. The Yangian Y(gl N ) is isomorphic to the tensor product Y(sIn) ® ZY(qI n ). 
See [1] for a proof. 

3.1.5. Let V be an irreducible finite-dimensional Y(g[^f)-module. Then there exists a unique 
vector v £ V such that 

T a b(u) v = 0, a > 6, 

T aa (V) v = Ca(u)v, a = l, ...,N, 

for suitable monic scries c a {v?j. IVIorcover, 

CaN P a (« + a) 



c a+ i(n) P a (u + a-l) 



1,...,N-1, (3.1) 



for certain monic polynomials P a {u). 

The polynomials Pi,...,Pjv-i are called i/ie Drinfeld polynomials of the module V. The 
vector t> is called a highest weight vector and the series ci(ii), . . . , cn(u) - the Yangian highest 
weights of the module V. 

For any collection of monic polynomials Pi , . . . , Pjv-i there exists an irreducible finite- 
dimensional Y(gl A r)-module V such that the polynomials Pi, ... , Pn-i are the Drinfeld polyno- 
mials of V. The module V is uniquely determined up to twisting by an automorphism of the 
form Xf- 

The claim follows from Drinfeld's description of irreducible finite-dimensional Y(slAr)-modu- 
les [2] and Proposition 13.11 

3.1.6. Let V%, V2 be irreducible finite-dimensional Y(gljv)-modules with respective highest 
weight vectors v\, V2- Then for the Y(gl7v)" m °dule V\ ® V2, we have 

T a b{u) v\ ® v 2 = 0, a > b, 

T aa {u) v x ®V2 = 4 ( n ) C ?H U ) vi®v 2 ., a=l,...,N. 

Let W be the irreducible subquotient of V\ ® V2 generated by the vector vi ® t^- Then the 
Drinfeld polynomials of the module W equal the products of the respective Drinfeld polynomials 
of the modules V\ and V2. 

3.1.7. A finite-dimensional irreducible Y(g[jv)-module V will be called polynomial if 

P N (u + N) 

cn{u) = p N ( u + N-iy (3 - 2) 

for some monic polynomial P/v(it). 
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For any collection of monic polynomials Pi, ... , P/v there exists a unique polynomial irre- 
ducible finite-dimensional ^(gt^-module V such that the polynomials P±, . . . , Pjy-i are the 
Drinfeld polynomials of V and (|3.2p holds. 

3.1.8. There is a one-parameter family of automorphisms 

p z : Y(gl N ) — ► Y(gl N ), T ab (u) h- T ab (u - z) , 

where in the right hand side, (u — z)~ l has to be expanded as a power series in u _1 . 

The Yangian Y(qI n ) contains the universal enveloping algebra U(qI n ) as a Hopf subalgebra. 
The embedding is given by the formula e a b i— > T 6a for all a, b. We identify [/(gljy) with its 
image. 

The evaluation homomorphism e : Y(qI n ) — > U(qI n ) is defined by the rule: T^ b t— > e ba for 

all a, 6 and T^ b 1 \—> for all a, 6 and all s > 1. 

For a gl^-module V denote by V(z) the Y (Ql N )-module induced from V by the homomor- 
phism e • p z . The module V(-z) is called the evaluation module with the evaluation point z. 

Let A = (A^ 1 ), . . . , A*-™-*) be a collection of integral dominant gt^-weights, where A® = 
(AP,...,A$) for t = l,...,n. For generic complex numbers . . . , 2^, the tensor product of 
evaluation modules 

M A (z) = M AC i) (Zl) ® • • • ® M ACn) (z„) 

is an irreducible Y(gljy)-module and the corresponding highest weight series c±(u), . . . , cn(u) 
have the form 

" „ _ z . + aW 

*m = n *_ + , • < 3 - 3 > 

i=l 4 
The corresponding Drinfeld polynomials are 

n A« 

^a(«)=n n («-^+s-a) 

for a = 1, . . . , N — 1. The Y(gljv)-module Ma. (2) is polynomial if A^ £ Z^o for ah Then the 
polynomial Pn(u) has the form 

A (i) 

1 = 1 s = l 

3.1.9. Consider as the gt^-module with highest weight (1, 0, . . . , 0). 

For any complex numbers z±, . . . ,z n , all irreducible subquotients of the Y(g[jv)-module 
C N (zi) ig> • • • <S> C N (z n ) are polynomial Y(g(^v)-modules. Moreover, for any polynomial irre- 
ducible finite-dimensional Y(g[7v)-module V, there exist complex numbers z\, . . . , z n such that 
V is isomorphic to a subquotient of the y(gl A r)-module C N (z±) <8> • • • <8> C N (z n ). 

The numbers z\, . . . ,z n are determined by the formula 

n N a-1 

H(u-z l )=niip a (u+s). 

i=l a=l s=0 

This formula follows from consideration of the action of the center of the Yangian Y(qI n ) in the 
module C^(zi) ® ■ ■ • ® C 7 ^). 
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3.1.10. A finite-dimensional 3^(fj^)-module will be called polynomial if it is the direct sum of 
tensor products of polynomial irreducible finite-dimensional ^(gl^-modules. 

If V is a polynomial finite-dimensional Y(0l7v)-module, then for any a,b = 1,...,N, the 
series T a b(u)\y converges to an End (F)-valued rational function in u. 

3.1.11. Let 7r : U(qI n ) — > End(C Ar ) be the representation homomorphism for the gl^y-modu- 
le C . Clearly, for any x E U(qI n ) we have 

[ir(x) <g> 1 + 1 ® x, T(u)) = 0. (3.4) 

For a non-degenerate matrix A € End(C^), define an automorphism va of Y(qI n ) by the 
formula 

(id®i/A)(T(u)) =J2 A~ l E ah A ® T ab (u). 

ab 

Let V be a finite-dimensional Yangian module with the representation [i : Y(qI n ) — > End (V) 
and /2 : GL^ — > End(V r ) the corresponding representation of the group GL^. The automor- 
phism induces a new Yangian module structure V A on the same vector space with the 
representation /i^ = /io z/4. Formula (|3.4p yields that for any x 6 y(gljy), 

/M(x) = £(A)/i(x)(/i(A))-\ (3.5) 

that is, the Y(g[jv)-modules V and V"" 4 are isomorphic. In particular, if V is a polynomial 
irreducible finite-dimensional Y(£|[jv)-module, then V A is a polynomial irreducible finite-dimen- 
sional Y(gl7v)- m odule too. 

3.2 Universal difference operator 

3.2.1. Let Q = (Qab) be an N x iV-matrix. Define 

N 

Xab(u,T~) = 5 ab QacTcbiu)^ 1 , a, b = 1, . . . , N. (3.6) 

c=l 

If V is a polynomial finite-dimensional Y(g( A r)-module, then X a b{u, r) acts on V"-valued functions 
in u, 

N 

f(u) l-> 5 ab f(u) - ^ QacT cb (u)f(u - 1). 
c=l 

Following [13], introduce the difference operator 

33(u,t)= ^ (-l) <T ^i CT1 («,r)-Y 2CT2 (u,r)---A' J v w («,T), (3.7) 

<tG5jv 

where the sum is over all permutations cr of {1, . . . , N}. The operator S)(it, r) will be called £/ie 
universal difference operator associated with the matrix Q. 

Lemma 3.2. Let ir be a map {1, . . . , N} — > {1, . . . , N}. If tt is a permutation of {1, . . . , N}, 

then 

J2 (- 1 )' 7 Xxinfar) X W2<72 {u,t) ■ ■ ■ X Wn<7n (u,t) = (-If £(u,r). 

(Tg5jv 

If it is not bijective, then 

^2 (-ly X w1(T1 (u,t) X W2(72 (u,t) ■ ■ ■ X wn(7n (u,t) = 0. 

The statement is Proposition 4.10 in [6]. 
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3.2.2. Introduce the coefficients To (it), . . . , T/v(it) of 2) (it, r): 

N 

3( V ) = E(-lf^(«)^ fc 

k=0 

The coefficients T&(ii) are called the transfer matrices of the XXX type model associated with Q. 
The main properties of the transfer matrices: 

(i) the transfer matrices commute: [Tfc(it), Tz(i>)] = for all k, I, u, v, 

(ii) if Q is a diagonal matrix, then the transfer matrices preserve the gljy- weight: [Tfc(it),e aa ]=0 
for all k, a, it, 

(hi) if Q is the identity matrix, then the transfer matrices commute with the subalgebra U(qI n ): 
[Tk(u),x] = for all k,u and x G U(gl N ), 

see [HIE]. 

3.2.3. Evidently, To (it) = 1. We also have Tat(u) = det Q qdet T(it), where 
qdetT(ix) = {-IT T lai {u)T 2(72 (u - 1) ■ ■ - T NaN (u - N + 1), 

and qdetr(u) = 1 + 0(it~ 1 ). 

Theorem 3.3. The coefficients of the series qdet T(u) are free generators of the center of the 
Yangian Y(gl N ). 

See [1] for a proof. 

3.2.4. If V is a polynomial finite-dimensional Y (gi N )-module, then the universal operator 
D(u,t) induces a difference operator acting on V- valued functions in u. This operator will be 
called the universal difference operator associated with Q and V and denoted by ®Qy(u, r). 
The linear operators %k(u)\v & End(y) will be called the transfer matrices associated with Q 
and V and denoted by %u Qv{ u )- They are rational functions in u. 

Example 3.4. Let 

V = M A (z) = M AC i) (zi) ® ■ ■ ■ ® M A(n) (z n ), 
Consider the algebra (^(st/v)) 18 " 1 - For a,b = 1, . . . ,N and z = 1, . . . , rt, define 

e« = ® e ab !•("-*) , L2(n, *) = 5 ah + -1- eg. 

The operator T {,(u) acts on Ma (2) as 

TV 

£ L<W~i(« - Zn) L^H^U ~ Zn-l) ■ ■ ■ L$ Cl {u - z 2 ) L« (u - Zl ). 

Cl,...,C„_l = l 

Lemma 3.5. If Dq(u,t) is the universal difference operator associated with the matrix Q and 
va : Y(q\ n ) — > Y(g[ N ) is the automorphism defined in Section 3.1.11, then 

va(Qq(u,t)) = Q aqa -i(u,t) (3.8) 

is the universal difference operator associated with the matrix AQA" 1 . 
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Proof. Consider matrices Q = (Q a b), T = (T a b(u)), and X = (X a b(u,T)). Then formula (|3.6p . 
defining X a b, may be read as X = 1 — QTt . Formula (|3.7j) for the universal difference operator 
may be understood as the row determinant of the matrix X , D(u, r) = det X. The definition of 
the automorphism ua reads as (id ® ua)(T) = A~ l TA. Then we have 

v a (Q(u,t)) = det (1 - QA~ 1 TAt- 1 ) = det - AQA^Tt^A). 

Now formula (|3,8p will be proved if we were able to write that last determinant as the product: 

det {A' 1 (l - AQA^Tt-^A) = det (A' 1 ) det (l - AQA^Tt^ 1 ) det A. 

The last formula may be proved the same way as the standard formula det MN = det M det N 
in ordinary linear algebra, using two observations. The first is that the entries of A are numbers 
and commute with the entries of X. The second observation is Lemma 13.21 describing the 
transformations of the row determinant of X with respect to row replacements. ■ 

3.2.5. Let V be a polynomial finite-dimensional Yangian module and fl : Gn — * GL(V) the 
associated GL^-representation. Then formulae (|3.5p and (|3.8p yield 

®AQA-i(u,t)\ V = KA)®q(u,t)\ v ^(A" 1 ). 



3.3 More properties of transfer matrices 

Let Q = diag (Qi, . . . , Qn) £ GL^. Let A = (A^, . . . , A( n )) be a collection of integral domi- 
nant gljv-weights, where AW = (A^\ . . . , A$) and A$ = for i = 1, . . . , n. 

For m = (mi, . . . , mjv) denote by M\(z)[m] C Ma(z) the weight-subspace of gl^r-weight m 
and by Sing M\(z)[m] C M\(z)[m] the subspace of g^-singular vectors. 

Note that the subspace Sing M\(z)[m] may be nonzero only if m is dominant integral, i.e. 
the integers mi, . . . , uin have to satisfy the inequalities m\ ^ • • • ^ rn/v- 

Consider the universal difference operator associated with Q and M\(z). Then the associated 
transfer matrices preserve M\(z)[m] and we may consider the universal difference operator 
®Q,M A (z)[m]( u i T ) acting on M\ (z) [m]- valued functions of u. We may write 

N 

^>Q,M A (z)[m](u,T) = ^2 ( _1 ) fc ^k,Q,M A ( Z )[m](u) T~ k . 
k=0 

As we know T ,Q,Af A (z)M (u) = 1, and we have 

^k,Q,MA{z)[m](u) = ^kO + ^klU~ l + 1fc2^~ 2 H 

for suitable T fci 6 End (Ma (z) [m]). 

Theorem 3.6. The fallowings statements hold. 

(i) The operators Tio, £20 > • • • > ^nq and Tn, T2I; ■ • • > ^Ni are scalar operators. Moreover, the 
following relations hold: 

N N 

x N +J2 (-1)" ^0 x N ~ k = H(x-Qi), 

k=l i=l 
N N N n 

k=l i=l 



x Qj 

3=1 
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(ii) Fork = l,...,N-l, 



N \ n ^-i , „ . aW 



^,q,m A (z)h(«) = n n 



u - z s + A)°> - i + 1 



u — 2 S — i + 1 

i=l / s=l j=l 

k n _, 



WaWmM =^fc( u )II n 



it - z s - i + 1 

1=1 s=l 

where Tfc(n) is a polynomial in u of degree nk. 

Proof. Part (i) follows from Proposition B.l in [6]. Part (ii) follows from the definition of the 
universal difference operator and the fact that the coefficients of the series T/v(u) belong to the 
center of Y(qI n ). ■ 

3.3.1. Assume that Q is the identity matrix. Then the associated transfer matrices preserve 
Sing M\(z)[m] and we may consider the universal difference operator '&Q=i j smgM A (z)[m]( u i t) 
acting on Sing M\ (z) [m]- valued functions of u. We may write 

N 

®Q=l,Sin g M A (z)lm](u,T)T N = £ (-l) k & k (u) (r - l) N ~ k (3.9) 

k=0 

for suitable coefficients &k(u). 

Note that the operators &k(u) coincide with the action in Sing M\(z)[m] of the modified 
transfer matrices Sk(u) from formula (10.4) of [6]. 

Theorem 3.7. The following three statements hold. 

(i) We have &q(u) = 1. 

(ii) For k = 1, . . . ,N, the coefficient &k(u) has the following Laurent expansion at u = oo: 

©k(u) = 6k,ou~ k + S^in-^- 1 + • • • , 

where the operators &i t o, ■ ■ ■ , &n,o are scalar operators. 
(Hi) For all d we have 

N N-k-1 N 

^(-l) fc 6 M J] (d-j) = H(d-m s -N + s). 

k=0 j=0 s=l 

Proof. Part (i) is evident. 

Since Q is the identity matrix, formula (|3.6p reads now as follows 

Xab(u) = 5 ab - T ab (u) T" 1 = (5 ab (r - 1) - uT 1 - 0{u~ 2 )) t~\ (3.10) 

Then part (ii) is straightforward from formulas (|3.9p and (|3.7p . 

Let v be any vector in Sing M\(z)[m] and d any number. To prove part (iii) we apply the 
difference operators in formula (|3,9p to the function vu d . The expansion at infinity of result of 
the application of the right side is 



N 



u d ~ N 



, k=0 



N-k-1 

n (d-j) 

3=0 



& kfi v + 0(u- 1 ) . 
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So it remains to show that 



2) 



Q=l,SingA/ A 



(z)[ m ](u,T)T N VU d = U d 



N 



~[(d - m s - N + s) 



+ o{u~ 1 ) 



as u goes to infinity. Since r N u d = u d (l + 0(u the last formula is equivalent to 



2> 



Q=l,Sing 



M A (z)[m]{u,T)vU d = U 



d-N 



N 



+ 0(u~ 1 )). (3.11) 



~[(d-m s -N + s) 

-S=l 

To prove ()3.1ip . observe that according to formula f|3. lOj) we have 

X ab (u) vu d = u d " 1 [[d5 ab - T^\u)) v + 0{u~ 1 )). 
Applying this remark to formula (|3.7p we get 

©0=1, Sing M A (z)[m] ( u > T ) v ^ 

£ ("I)' 7 {(d-N+ 1) <5 Vl - lg> ) • • • {dS N>(TN - r«J v + O^ 1 ) ] . 



d-N 



Each element T^ b acts as ^ ej^ in M\{z), which corresponds to the standard gl^-action in 

i=i 

the g^-module M\. Since i> is a {^-singular vector of a gl^r-weight (mi, . . . , mjv) ; we have 

v 



X; (-i) CT ((rf - jv + 1) s hai - jfa • • • (ds N , aN - 



~[(d - m s - N + s) 



v. 



Indeed only the identity permutation contributes nontrivially to the sum in the left side. This 
proves part (iii) of Theorem 13.71 ■ 

3.4 First main result 

Theorem 3.8. Let V be a polynomial finite- dimensional Y(gl N )-module and Q £ GLkt. Con- 
sider the universal difference operator S)q,v(u,t) associated with Q and V . Then the kernel of 
®Q t y(u,T) is generated by quasi- exponentials. 

The theorem will be proved in Section [5j 



4 Continuity principle for difference operators 
with quasi- exponential kernel 

4.1 Independent quasi-exponentials 

Let V be a complex vector space of dimension d. Let p £ C[u] be a monic polynomial of degree k. 
Consider the differential equation 



for a V- valued function f(u). Denote by W p the complex vector space of its solutions. The map 

S P : W p - V®\ w i ^ {w(0),w'(0),...,w^- l \0)), 
assigning to a solution its initial condition at u = 0, is an isomorphism. 
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Let Ai, . . . , A; be all distinct roots of the polynomial p of multiplicities k\, . . . , ki, respectively. 
Let vi, . . . , Vd be a basis of V . Then the quasi-exponentials 

e x i u u a v b , j = l,...,l, a = 0,...,kj-l, b = l,...,d, (4.1) 

form a basis in W p . 

Lemma 4.1. Assume that X a — Xb 2-7riZ for all a ^ b. Then the kal functions listed in (|4.ip 
are linear independent over the field of one-periodic functions. 

4.2 Admissible difference operators 

Let Aq(u), . . . , Aj\f(u) be End (y)-valued rational functions in u. Assume that each of these 
functions has limit as u — > oo and Aq(u) = 1 in End(y). For every k, let 

A k (u) = A% + A^u- 1 + A% 2 u- 2 + ■■■ 
be the Laurent expansion at infinity. Consider the algebraic equation 

det (^tv,o + ^-1,0 + • • • + x N ~ 1 A lfi + x N ) = (4.2) 
with respect to variable x and the difference operator 

N 

V = ^2A k (u)T- k 

k=0 

acting on ^-valued functions in u. Equation (|4.2j) will be called the characteristic equation for 
the difference operator 2). 

The operator T> will be called admissible at infinity if det A/v,o 7^ 0) or equivalently, if x = 
is not a root of the characteristic equation. 

Example 4.2. Let V be the tensor product of polynomial finite-dimensional irreducible Y(qI n )- 
modules and Q € GL^. Let ®Qy(u, r) be the associated universal difference operator. Then 
®q,v(u,t) is admissible at infinity, see Sections 3.2.3 and !3.3[ 

Lemma 4.3. Assume that D is admissible at infinity and a nonzero V -valued quasi- exponential 
Q u (u d Vd + u d ~ l Vd-i + • • • + vq) lies in the kernel of 2). Then Q is a root of the characteristic 
equation (|4.2D . 

4.3 Continuity principle 

Let Ao(u,e), . . . ,A^(u,e) be End (V)-valued rational functions in u analytically depending on 
e € [0, 1). Assume that 

TV 

• for every e € [0, 1) the difference operator D e = ^ A k (u, e)r~ k is admissible at infinity, 

k=0 

• for every e € (0, 1) the kernel of D e is generated by quasi-exponentials, 

• there exists a natural number m such that for every e € (0, 1) all quasi-exponentials 
generating the kernel of D t are of degree less than m. 

Theorem 4.4. Under these conditions the kernel of the difference operator 2) e= o is generated 
by quasi-exponentials. 
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Proof. For every e £ [0,1), the characteristic equation for T> e has Nd roots counted with 
multiplicities. As e tends to the roots of the characteristic equation of D e tend to the roots of 
the characteristic equation of %D t =o- All these roots are nonzero numbers. For small positive e 
the set of multiplicities of roots does not depend on e. 

The following lemma is evident. 

Lemma 4.5. There exist 

• a number e with < e ^ 1, 

• for any e, < e < e, a way to order the roots of the characteristic equation ofT) t (we 
denote the ordered roots by Q\, . . . , Q%d)j 

• a way to assign the logarithm to every root Q e j 

such that for every j the number g| continuously depends on e and = qf whenever Q e - = Q\- 

Let m be the number described in Section [4.31 For every e, < e < e, we define p e € C[u] to 
be the monic polynomial of degree k = mNd, whose set of roots consists of m copies of each of 
the numbers q\, . . . , q e Nd . 

Let W Pe be the fed-dimensional vector space of quasi-exponentials assigned to the polyno- 
mial p e in Section B~T1 By assumptions of Theorem 14.11 for every e, < e < e, the space W Pt 
contains an iV(i-dimensional subspace U e generating the kernel of D e . This subspace determines 
a point in the Grassmannian Gr(W Pt , Nd) of A^d-dimensional subspaces of W Pe . 

The map p Vl : W Pl — > V® k identifies the Grassmannian Gr(W Pt , Nd) with the Grassmannian 
Gi(V® k , Nd). The points p Pc (U e ) all lie in the compact manifold Gi(V® h ,Nd) and the set 
of all such points has an accumulation point U € Gr(V® k , Nd) as e tends to zero. Then 
P P } =0 (U) C Gr(Wp e=0 , Nd) is an iV<i-dimensional subspace of V- valued quasi-exponentials. Using 
Lemma l4~T| we conclude that the space p p ^ =0 (U) generates the kernel of S) e =o- B 

5 Bet he ansatz 
5.1 Preliminaries 

Consider as the gl^-module with highest weight (1, 0, . . . , 0). For complex numbers z\, . . . , 
gz n , denote z = (z\, . . . , z n ), and 

M{z) = C N ( Zl )® ■■■®C N (z n ), 

which is a polynomial Y(gi N ) -module. Let 

M(z) = ® mi ,..., mN M(z)[mi, . . .,m N ] 

be its q [^-weight decomposition with respect to the Cartan subalgebra of diagonal matrices. 
The weight subspace M(z)[mi, . . . , rajy] is nonzero if and only if mi, . . . , mjv are non-negative 
integers. 

Assume that Q = diag (Qi, ■ . . , Qn) is a diagonal non-degenerate iV x iV-matrix with distinct 
diagonal entries, and consider the universal difference operator 

N 

®q,m(z)(u,t) = ^(-l) fc T fciQiM(z) (n)r- fc . 

k=0 

associated with M(z) and Q. Acting on M(z)-valued functions the operator i S)q : m(z)( u i t ) 
preserves the weight decomposition. 

In this section we shall study the kernel of this difference operator, restricted to M(z)[mi, . . . , 
m7v]-valued functions. 
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5.2 Bethe ansatz equations associated with a weight subspace 

Consider a nonzero weight subspace M(z)[m\, . . . , mjy]. Introduce I = . . . , Zjv-i) with lj = 
nij + i + • • • + rriM- We have n ^ l\ ^ • • • ^ Zjv— l ^ 0. Set Iq = In = and / = ?! + ••• + Zjv_i. 
We shall consider functions of I variables 

f _ (JXi f W f (2) ,(2) .(JV-1) .(iV-l)N 

i — j • • • , % j <a j • • • i % , • • • j H j • • • j tj JV _ 1 j • 

The following system of / algebraic equations with respect to I variables t is called the Bethe 
ansatz equations associated with M(z)[mi, . . . , m^r] and Q, 

n h h 

qi n - *. + 1) n - ^ - 1) n - f ) ^ 

s=l j'=l j'=l 

n h h 

s=l j'=l j'=l 

Qa n (tf - tf- x) + 1) n {tf - tf ~ i) n (*« - 

j'=i ^'=1 j'=i 

= Q a+1 n {tf - tf- i} ) n (t } - # } + o ff d a) - - i ) > 

i'=i j'=i i'=i 

In-2 In-i 

qn-i n (4 a) - + !) n (4 a) - f ? - x ) 
j'=i f=i 

In-2 In-i 

(a) ,(a-l)\ TT /.(a) .(a) 



Qtv] [ (tf-tf~ l) ) ] [ (tf-tf + 1). 



j'=l j'=l 

Here the equations of the first group are labeled by j = 1, . . . , Zi, the equations of the second 
group are labeled by a = 2, . . . , N — 2, j = 1, . . . , l a , the equations of the third group are labeled 
by j = 1, .. ■ ,Zjv-i- 

A solution i of system (|5.ip will be called off-diagonal if f^- 7^ for any a = 1, . . . , N — 1, 
1 < j < j' < ^, and / for any a = 1, . . . , N - 2, j = 1, . . . , l a , f = 1, . . . , l a+1 . 

5.3 Weight function and Bethe ansatz theorem 

Denote by u(t, z) the universal weight function associated with the weight subspace M(z)[mi, 
. ..,mjv]- The universal weight function is defined in formula (6.2) in [6], see explicit for- 
mula (15.21) below. At this moment, it is enough for us to know that this function is an 
M(z)[mi, . . . , mAr]-valued polynomial in t, z. 

If t is an off-diagonal solution of the Bethe ansatz equations, then the vector uj(t, z) € 
M(z)[mi, . . . , mjv] is called the Bethe vector associated with t. 

Theorem 5.1. Let Q be a diagonal matrix and t an off-diagonal solution of the Bethe ansatz 
equations (15. lj) . Assume that the Bethe vector u(t, z) is nonzero. Then the Bethe vector is an 
eigenvector of all transfer-matrices ^k,Q,M(z)( u ); k = 0,. . . ,N. 

The statement follows from Theorem 6.1 in [6]. For k = 1, the result is established in [3]. 
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The eigenvalues of the Bethe vector are given by the following construction. Set 



+ lA«-f"l 



x\u, t, z)=q 1 h n —^i 

j=1 U-tj +1 j=1 U-tj 
for a = 2, . . . , N. Define the functions Xk(u, t, z) by the formula 

N 

(1 - X \u, t, z) r" 1 ) • • • (1 - X N (u, t, z) r" 1 ) = £ (-l) fc X k (u, t, z) r~ k . 

Then 

= X k (u,i,z)ui{i,z) 
for k = 0, . . . ,N, see Theorem 6.1 in [6]. 

5.4 Difference operator associated with an off-diagonal solution 

Let t be an off-diagonal solution of the Bethe ansatz equations. The scalar difference operator 

JV 



k=0 

will be called the associated fundamental difference operator. 

Theorem 5.2. The kernel ofDz(u,T) is generated by quasi- exponentials of degree bounded from 
above by a function in n and N . 

This is Proposition 7.6 in j9], which is a generalization of Proposition 4.8 in [8]. 
5.5 Completeness of the Bethe ansatz 

Theorem 5.3. Let z\, . . . ,z n and Q = diag (Qi, ■ ■ ■ ,Qn) be generic. Then the Bethe vectors 
form a basis in M(z)[m\, . . . ,m^]. 



Theorem 15.31 will be proved in Section 15.61 
Corollary 5.4. Theorems 15.31 and 14.41 imply Theorem 13.81 

Proof. Theorems 15.21 and 15.31 imply that the statement of Theorem 13.81 holds if the tensor prod- 
uct M{z) is considered for generic z and generic diagonal Q. Then according to the remark 
in Section 3.2.5, the statement of Theorem 13.81 holds if the tensor product M{z) is considered 
for generic z and generic (not necessarily diagonal) Q. Then the remark in Section 3.1.9 and 
Theorem 14.41 imply that the statement of Theorem 13.81 holds for the tensor product of any 
polynomial finite-dimensional irreducible y(g[jy)-modules and any Q € GL^. Hence the state- 
ment of Theorem 13.81 holds for direct sums of tensor products of polynomial finite-dimensional 
irreducible y(g [^-modules and any Q € GLjy. ■ 
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5.6 Proof of Theorem 15.31 

5.6.1. For a nonzero weight subspace M(z)[mi, . . . , mjv] C M(z) denote by d[m\, . . . , mjv] its 
dimension. Let n ^ l\ ^ • • • ^ ZjV-i ^ be the numbers defined in Section I5T21 

A vector a = (ai, . . . ,a n ) with coordinates aj from the set {0, 2, 3, . . . , iV} will be called 
admissible if for any j = 1, . . . , N — 1 we have lj = Gtj \ i = 1, . . . , n, and a, > J } In other 
words, a is admissible if rrij = a, | i = 1, . . . , n, and a, = _?'}. 

If a = (oi, . . . , a n ) is admissible, then for j = 1, . . . , N — 1, there exists a unique increasing 
map p a>i : {1, . . . , k} -> {1, . . . , n} such that #p~}(j) = 1 if % > z and #/9~)(j) = if a, < i. 

We order admissible vectors lexicographically: we say that a > a' if a^v = ctjy, ajv-i = a N-i' 
... , ai = a[, ai-x > ai-\ for some i. 

5.6.2. Let v = (1,0, ... ,0) £ be the highest weight vector. Consider the set of vectors 
e a v = e ai) iv ® ■ • • ® £a n ,iv £ M(z) labeled by admissible indices a. Here e au \v denotes v if 
ai = 0. Then this set of vectors is a basis of the weight subspace M{z)\m\, . . . , mjv]- In 
particular, the total number of admissible indices equals d[mi, . . . , m^]- 

5.6.3. For a function f(u\, . . . , ti^) set 
^y m ni,...,M fe 

...,%)= ^ /(n CT1 , . . . ,u CTfe ), 

where the sum is over all permutations a of {1, . . . , k}. 

Lemma 5.5. The universal weight function oj{t,z) is given by the rule 



n h N-l h-i h 

:«.«) = n n cf -*) n nncf-f) 

s=l j=l 6=2 i=l j=l 



(5.2) 



(i) /i) ■ ■ • Sym (jv-i) (jv-i) 



n / n s-1 — Z r + L 



s=l 
a 3 >l 



(1) _ 11 (1) _ 



JV-l n 

n n 



6=2 



1 \ l „-i r„\ ~ r „-i 



n 



Pa,b( S > Pa,b~l( r ) 



Ab) (6-1) 
r=l l „-l^i t „-l 



where the sum is over all admissible a . 



The lemma follows from formula (|3.3p . formula (6.2) in [B], and Corollaries 3.5, 3.7 in |14j . 

5.6.4. Let £i, . . . , z n be real numbers such that Zi+\ — z% > iV. Assume that Q — diag (1, q, . . . , 
q N ~ l ), where q is a nonzero parameter. Then the Bethe ansatz equations (|5,ip take the form 



IK* 



m-,, + i ) n (t m- f <:)- 1) n(4 1> -4?) 

i'=i i'=i 



(5.3) 



i'=i 
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n {tf ] - t<r i} + 1) n {tf - tf - 1) n ^ - # +i >) 
j>=i f=i f=i 

j'=i j'=i j'=i 

^iV-2 'jV-1 'jV-1 

n (t } - + !) n (4 a) - - j ) = « n (*f - n (4 a) - ^ + 
j'=i /=i j'=i /=i 

For g = the right hand sides of equations (|5.3|) equal zero and to solve the equations for q = 
one needs to find common zeros of the left hand sides. 

Let a = (ai, . . . , ajv) be an admissible index. Define t(a,0) = (t^'\a,0)) to be the point in 

C' with coordinates $p(a,, 0) = z Pa .q) — i for all i, j. Then t(a, 0) is an off-diagonal solution of 
system (|5.3p for g = 0. That solution has multiplicity one. 

Hence for every small nonzero q, there exists a unique point i(a,q) = (t^\a,q)), such that 

• t(a, q) is an off-diagonal solution of system (j5.3f) with the same q, 

• t(a,q) holomorphically depends on q and tends to t(a, 0) as g tends to zero. 

Therefore, for all i, j, we have 

^ ) (a,g) = ^ (i) -i + 0(g). (5.4) 
Lemma 5.6. As q tends to zero, the Bethe vector u(t(a,q), z) has the following asymptotics: 

u(t(a,q),z) = C a e a v + 0{q) H , 

where C a is a nonzero number and the dots denote a linear combination of basis vectors e a >v 
with indices a' lexicographically greater than a. 



The lemma follows from formulae (|5.4p and (|5.2() , 
Lemma 15.61 implies Theorem 15.31 

6 Comparison theorem 

6.1 The automorphism Xf an d the universal difference operator 

Let f(u) be a rational function in u whose Laurent expansion at u = co has the form f(u) = 
1 + 0(n _1 ). Then the map x.f '■ T(u) h- > f(u)T{u) defines an automorphism of the Yangian 
Y(qI n ), see Section 3.1.4. 

For a Yangian module V, denote by V* the representation of the Yangian on the same vector 
space given by the rule 

X\yf = iXf{X))\ V 

for any X G Y(gl N ). 

Lemma 6.1. Let V be an irreducible finite- dimensional Y(Ql N )-module with highest weight 
series c\(u), . . . ,cn(u). Then the Yangian module V? is irreducible with highest weight series 
/(u)ci («),... ,f(u)c N (u). 
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Lemma 6.2. Let V be a finite- dimensional polynomial Y(gl N ) -module. Then the corresponding 
transfer matrices satisfy the relation 

Za(u)\ V f = f(u) ■■■/(«- a + 1) % a (u)\ v 
fora = l,...,N. 



6.2 Comparison of kernels 

Let V be a finite-dimensional polynomial y(g[jv)-module. Consider the universal difference 
operators ®Qy(u, r) and Dq yf(u, r). Let C(n) be a function satisfying the equation 

C(u) = f{u)C(u-l). 

Then a ^-valued function g(u) belongs to the kernel of DQy(u,r) if and only if the function 
C(u)g(u) belongs to the kernel of Dq v j{u,t). 

6.2.1. Let A = (Ai,...,Ajv) be a dominant integral gt^-weight and a a complex number. 
Denote A(a) = (Ai + a, . . . , Ajy + a). Let V = M\{z) be an evaluation module and 

fw 



u — z — a 



Then V? = M A{a) (z + a), see Section 3.1.8. 

Let A = (A^ 1 ), . . . , A^" - )) be a collection of integral dominant gl^-weights, where A^ = 
(Aj ,...,A^') for % = 1, ...,n. Let a = (ai,...,a n ) be a collection on complex numbers. 
Introduce the new collection A(o) = (A( 1 )(a 1 ), . . . , A^ n \a n )). 

For complex numbers z±, . . . , z n , consider the tensor products of evaluation modules 

M A (z) = M A( i) (zi) <g> • • • ® M A(n) (z„), 

M A(o )(z + a) = M A( i )(ai) (zi + ai) ® • • • <g> M A(n ) (0n) (z n + a n ). 
Let F = M A (z) and 

/(«) = n 



u - z s 



U-Z s 
s=l 



Then V f = M A(a) (z + a) . 

Theorem 6.3. An M \(z) -valued function g(u) belongs to the kernel of ?Oq 1 m a (z)(. u i t ) an d 
only if the function C(u)g(u) belongs to the kernel of Dq^,j ^ z+ci ^(u,t) , where 



c ^ = n f 



T(u -zs + l) 



L T(u - z s - a s + 1) ' 

6.2.2. Assume that M A (z) is a polynomial Yangian module and a = (ai,...,a/v) are non- 
negative integers. Then C(u) is a polynomial and M A ( a )(2: + a) is a polynomial Yangian module. 
In this case Theorem 16.31 allows us to compare the quasi-exponential kernels of the difference 
operators D Q Ma{z) (u, t) and D Q:MA{a)[z+a) (u,T). 
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7 The kernel in the tensor product of evaluation modules 

7.1 The second main result 

Let Q = diag (Q\, . . . ,Qn) £ GLn- Let A = (A^ 1 ), . . . , A^) be a collection of integral dominant 
gl^-weights, where AW = (Aj 1 , . . . , A^j) and A^ = for i = 1, . . . , n. 

Choose a weight subspace M/y(z)[mi, . . . , m^] C M\(z) and consider the universal difference 
operator ^>Q.M A (z)[m] associated with Q and M\(z)[mi, . . . , mjv]- 

Theorem 7.1. Assume that the nonzero numbers Q±, ■ ■ ■ ,Qn are distinct and the argument of 
each of them is chosen. Then for any i = 1, . . . , N and any nonzero vector vq € M\(z)[mi, . . . , 
mjv]j there exists an M\{z)\m\, . . . ,m]\f]-valued quasi- exponential 

Q?(v u m > + Vl u^- 1 + ■■■ + v mi ) (7.1) 

which lies in the kernel ofDQ t M A (z)\m]- Moreover, all such quasi- exponentials generate the kernel 

0f2>Q,M A (z)[m]- 

Proof. On one hand, by Theorem l3.8l the kernel of 2Q : M A (2)[m] is generated by quasi-exponen- 
tials. On the other hand, by Lemma 14.31 and part (i) of Theorem 13.61 an Y quasi-exponential 
lying in the kernel of 3q,m a (2)H must be of the form (|7.1|) . where vq is a nonzero vector. Since 
the kernel is of dimension N ■ dim M\(z)[m], there exists a quasi-exponential of the form (|7.ip 
with an arbitrary nonzero Vq. ■ 

7.1.1. For i = 1, . . . , n, denote Si = {z { - 1, - A$_i - 2, . . . , z { - A® - N}. 

Theorem 7.2. Assume that index i is such that Z{ — zj £ Z for any j ^ i. Let f(u) be an 
M\(z)[m]-valued quasi- exponential lying in the kernel of ^Q t M A (z)[m\- Then 

(i) f(u) is uniquely determined by its values 

v N = f(zi - 1), v N -! = f(zi - A$_ a - 2), . . . , vi = f{zi - A® - N) e M A (z)[m} 
at the points of Si . 

(ii) If vn = vn-i = ■ ■ ■ = v j = for some j > 1, then f(z{ — k) = for k = 1,2,..., 
A^ + iV-i + 1. 

(in) For any vectors vn,vn-i, ■ ■ ■ , v i £ M\(z)[m], there exists a quasi- exponential f(u) which 
lies in the kernel of T>Q t M A (z)[m] an d takes these values at Si- 

Proof. Consider the polynomial difference operator 

n N-l 

S = II Y[( u - z s~ J + 1 ) ®Q,M A (z){m] 
n N-l n N-l 

= II - Zs ~ 3 + 1) + ~[[~[[(u - z s - j + lj'Z^T- 1 + ■ ■ ■ 

s=l 3=1 s=l j=1 



+ Y[(u-z s -N + 2) i N - 2 (u) t~ n+2 + %n-i{u) t- 

s=l 

N n N-l 

+ (-if n^niK" - A S s) - 1 + i ) t ~ n > 

i=l s=l i=l 

see notation in part (i) of Theorem 13.61 We have (T)f)(u) = 0. 



20 



E. Mukhin, V. Tarasov and A. Varchenko 



Condition (Df)(zi + N — 2) = gives the equation f(zi — 2) = const f(zi — 1) and thus the 
value f(zi — 1) determines the value f(zi — 2). Condition (T) f)(zi + N — 3) = gives the equation 
/(zj — 3) = const f(zi — 2) and hence the value f(zi — 2) determines the values f(zi — 3). We 
may continue on this reasoning up to equation f(zi — A^_ x — 1) = const f(zi — A^_ x ) which 
shows that the value f(zi — A®^) determines the value f(zi — A^_ 1 — 1). 

Condition 0bf){zi + N — A^_ x — 2) = does not determine the value f(zi — A^_ 1 — 2), since 
f(zi~ Ajy^— 2) enters that condition with coefficient 0. But condition 0bf)(zi+N— A$i— 3)=0 
gives the equation /(zj — A^'_ 1 — 3) = const /(zj — A^_j — 2) which shows that the value 
f(zi — A^_ x — 3) is determined by the value f(z-i — — 2). 

Now we may continue on this reasoning up to the equation f(zi — A^_ 2 — 2) = const f{zi — 

A^_ 2 — 1) which determines f(zi — A$_ 2 — 2) if f{zi — A^_ 2 — 1) is known. Condition + 

N — A$_ 2 — 3) = does not determine the value f(z{ — A®_ 2 — 3), but condition {J)f)(z{ + 

N — A^_ 2 — 4) = gives the equation f(zi — A^_ 2 — 4) = const f{zi — A^_ 2 — 3) and so on. 
Repeating this reasoning we prove parts (i) and (ii) of the theorem. 

The same reasoning shows that if f(u) = for u £ Si, then the quasi-exponential f(u) 
identically equals zero. Since the kernel of D is generated by quasi-exponentials, we obtain part 
(hi) of the theorem. ■ 

7.1.2. Assume that Q is the identity matrix. Consider the subspace Sing M\(z) [m] C M\(z) [m] 
of gljy-singular vectors and the associated universal difference operator 53q =1j g ing M A (z)[m] ( u > T )- 

Theorem 7.3. For any i = 1, . . . , N and any nonzero vector vo G Sing M\(z)[mi, . . . ,m^], 
there exists a Sing M\(z)[mi, . . . ,m]\f]-valued polynomial 

VQ u m i+ N-i + vi u m l+ N-i-l + . . . + Wjv _. (7 . 2 ) 

which lies in the kernel ofT>Q = i^- mg M A {z)\m] • Moreover, all such polynomials generate the kernel 

Of £>Q=1, Sing M A {z)[m] ■ 

Proof. On one hand, by Theorem 13.8} the kernel of 3)q=i, sing M A (z)[m] is generated by quasi- 
exponentials. On the other hand, by part (i) of Theorem 13.61 any quasi-exponential lying 
in the kernel of 3)Q=i ) singAfA(*)[m] m ust a polynomial. By Theorem 13.71 such a polynomial 
has to be of the form indicated in (|7.2p . where vq is a nonzero vector. Since the kernel is of 
dimension iV-dim Sing M\(z)[m], there exists a polynomial of the form (|7.2p with an arbitrary 
nonzero vq. ■ 

7.2 Kernel of the fundamental difference operator associated 
with an eigenvector 

Assume that v € M\(z)[m] is an eigenvector of all transfer matrices, 

^Q,M A (z)[ m ](u)v = \,v(u)v, k = 0,... ,N. 

Then the scalar difference operator 

N 



V v (u,t u ) =^(-l) fc A M (u)r- 



k=0 



will be called the fundamental difference operator associated with the eigenvector v. 

Theorems 17.11 17.21 and 17.31 give us information on the kernel of the fundamental difference 
operator. 
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Corollary 7.4. 

(i) Assume that Q G GLn is diagonal with distinct diagonal entries. Then the kernel of 
^>v(u, t u ) is generated by quasi- exponentials Qfpi(u), . . . , Q n Pn{u), where for every i the 
polynomial pi(u) G C[u] is of degree mi. 

(ii) Assume that Q is the identity matrix and the eigenvector v belongs to Sing M\(z)[m]. 
Then the kernel of® v (u, t u ) is generated by suitable polynomials pi(u), . . . ,pn(u) of degree 
mi + N — 1, . . . , mjv, respectively. 

Corollary 7.5. Assume the index i is such that Zi — Zj ^ Z for any j ^ i. Let f(u) be a quasi- 
exponential lying in the kernel of D v (u,t u ). Then 

(i) f(u) is uniquely determined by its values 

v N = f( Zi -l), v N - 1 = f(z i -A$_ 1 -2), v l = fizi-kf -N) 

at the points of Si . 

(i) 

(ii) If vn = vn-i = ■ • • = v j = for some j > 1, then f{Zi — k) = for k = 1, 2, . . . ,Aj_j_ + 
N-j + 1. 

(Hi) For any numbers vn, vn-i, ■ ■ ■ , v\ there exists a quasi- exponential f(u) which lies in the 
kernel of D v (u,t u ) and takes these values at Si. 

If v is a Bethe eigenvector, then these two corollaries were proved in [8] and [9]. 

8 Generating operator of the Gaudin transfer matrices 
8.1 Current algebra flljvf 33 ] 

8.1.1. Let 0[jv[a?] be the Lie algebra of polynomials with coefficients in gi N with point-wise 
commutator. The elements = e a bx s with a, b = 1, . . . , N, s = 0, 1, . . . , span gl^-fx]. We 

have [e^,e^] = o~bc e ^ S ^ — ^ad^cb +S ^ ■ We shall identify qI n with the subalgebra of flljvM of 
constant polynomials. 

For a, b = 1, . . . , N, we set 

oo 

m«) = !>£?«— 1 

and L{u) = E a ,b E ab®L ab {u) G End (C*) <8> gl^K^" *))• 

Let V be a gt^-module. For z G C denote by M(z) the corresponding evaluation gljvM - 
module on the same vector space, where we define e^ b \v = z s e a b for all s, a, b. Then the series 
Lab(u)\v converges to the End (y)-valued function (u — z)~ l e ba , which is a rational function 
in u. 

Let A = (Ai, . . . , Ajv) be an integral dominant gl^-weight, M\ the corresponding irreducible 
highest weight gl^r-module, Ma(z) the associated evaluation qI n [x]-module. Then the module 
M\(z) is called polynomial if An is a non-negative integer. 

8.1.2. Let A = (A^ 1 ), . . . , A^ n )) be a collection of integral dominant gl^-weights, where A^ = 
(Aj*\...,A^) for i = l,...,n. For zi,...,z n G C, we consider the tensor product of gl^fx]- 
evaluation modules: 

Ma(z) = M A(1) (zi) <g> ■ •• ® M Ai n)(z n ). 
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For any a, b, the series L ab (u) acts on M\(z) by the formula 

n 

Lab{u)\ MA (z) = Yl 

j=l 3 

If zi, . . . , zn are distinct, then A4a(z) is an irreducible gl^xj-module. Let 

M AC i) <8> • • • <8> M A ( n ) = 0a M A 

be the decomposition of the tensor product of gl^-modules into the direct sum of irreducible 
gljy-modules. Then for any z € C, 

M A( i) (*) <g> • • • <g> M A( „) (z) = © A M A (z) 

is the decomposition of the tensor product of evaluation gl^fx] -modules into the direct sum of 
irreducible g^fx] -modules. 

8.1.3. Consider as the gt^-module with highest weight (1, 0, . . . , 0). Let A = (Ai, . . . , An) 
be an integral dominant gl^r-weight with A^r € Z^o- Then there exists k £ Z^o such that (£, N )® k 
contains M A ELS cl gljy-submodule. 

The previous remarks show that for any Zi,...,z n E C, all irreducible submodules of the 
gl^x] -module C N (zi) (g) ■ ■ ■ <g> C N (z n ) are tensor products of polynomial evaluation gl^-fx]- 
modules. Moreover, if V is a gl^rfx] -module which is the tensor product of polynomial eval- 
uation g^xj-modules, then there exist z\, . . . , z n such that V is isomorphic to a submodule of 
the gl N [x] -module C N (z{) <8> ■ • • ® C N (z n ). 

8.1.4. Let 7r : U(qIn) — > End (C^) be the representation homomorphism for the gl^r-module C^. 
Clearly, for any x £ U(gl N ) we have 

[tt(x) <g> 1 + 1 <g> x, L(it)] = 0. (8.1) 

For a non-degenerate matrix ^4 € End(C Ar ), define an automorphism v& of g[j\r[x] by the 
formula 

(id®v A )(L(u)) = A- l E ab A®L ab (u). 

ab 

Let V be a finite-dimensional flljvM -module with the representation fi : gljvfx] — » End (V) 
and /i : GL^v — > End (V) the corresponding representation of the group GL^. The auto- 
morphism induces a new gljy [x] -module structure V A on the same vector space with the 
representation [la = H° va- Formula (18. ip yields that for any x E gl N [x], 

that is, the g [x] -modules V and V A are isomorphic. In particular, if V is the tensor product of 
polynomial evaluation g[^y[x] -modules, then V A is the tensor product of polynomial evaluation 
0ljv[x] -modules too. 

8.2 Fundamental differential operator 

Let K = (K ab ) be an N x N matrix with complex entries. For a, b = 1,...,N, define the 
differential operator 

X a b{u, d u ) = 5 ab d u - K ab - L ab (u), 

where d u = d/du. 
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Following [13], introduce the differential operator 

V(u,d u )= ^2 ( _1 ) CT x iai(u,d u ) X 2 a 2 (u,d u ) ■ ■ ■ X NaN (u,d u ), 

where the sum is over all permutations a of {1, . . . , N}. The operator T>(u,d u ) will be called 
the universal differential operator associated with the matrix K. 

Lemma 8.1. Let tt be a map {1, . . . , N} — > {1, . . . , N}. If tt is a permutation of {1, . . . , N}, 

then 

Yl X 7T 1 a l (u,r) X n2U2 (u,T) ■ ■ ■ X nNaN (u,T) = (-l) n V(u,d u ). 

If tt is not bijective, then 

Yl ( -1 ) CT X tio-i( u ' t ) X 7T2o- 2 (%t) ■ ■ ■ X 7r N a N (u,r) =0. 

The statement is Proposition 8.1 in [6]. 

8.2.1. Introduce the coefficients So(u), . . . ,Sn(u) of T>(u, d u ): 

N 

V(u,d u )=Y J (-V k S k {u)d^-\ 

k=0 

in particular, So(u) = 1. The coefficients Sk(u) are called the transfer matrices of the Gaudin 
type model associated with K. 

The main properties of the transfer matrices: 

(i) the transfer matrices commute: [Si~(u),Si(v)] = for all k, I, u, v, 

(ii) if K is a diagonal matrix, then the transfer matrices preserve the q [^-weight: [Sk(u), e aa ]=0 
for all k, a, u, 

(iii) if K is the zero matrix, then the transfer matrices commute with the subalgebra U{q\^): 
[Sk(u),x] = for all k, u and x G U(gl N ), 

see [HIE]. 

8.2.2. If V is the tensor product of evaluation finite-dimensional 3 t^[x] -modules, then the 
universal operator T>(u, d u ) induces a differential operator acting on V-valued functions in u. 
This operator will be called the universal differential operator associated with K and V and 
denoted by T>K,v{u,d u ). The linear operators Sk(u)\v € End (V) will be called the transfer 
matrices associated with K and V and denoted by k vi u )- They are rational functions in u. 

8.2.3. If T>k(u,t) is the universal differential operator associated with the matrix K and 
is A : [jv [ x ] ~~ * 0'v[ x ] i s the automorphism defined in Section 8.1.4. Then Lemma 18.11 implies 
that 

iy A {1^K(u,d u )) = T>AKA-i(u,d u ) 

is the universal differential operator associated with the matrix AKA" 1 , cf. Lemma 13.51 

Let V be the tensor product of finite-dimensional evaluation gl^fzj-modules, and \x : Gn — > 
GL{V) the associated GL^-representation. Then 

V AKA-A u ,du)\v = fi{A)V K (u,d u )\ v fi(A~ l ). 
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8.3 More properties of the Gaudin type transfer matrices 

Let K = diag (Ki, ... ,Kj\f). Let A = (A^ 1 ', . . . , A^ n ^) be a collection of integral dominant 
gljv-weights, where = (A^ , . . . , A$) and A$ = for i = 1, . . . , n. 

For m = (mi,...,mjv) denote by A4\_{z)[m] C M.\(z) the weight subspace of the fly- 
weight m and by Sing M\{z)\m\ C M.A{z)[m\ the subspace of gt^-singular vectors. 

Consider the universal differential operator 

TV 

Vk,m a (z)(u,9u) = ^(-l) k Sk,K,M A {z)(u)d^- k 

k=0 

associated with K and M\(z). We have So,k,Ma(z)\™](. u ) = an d 

<Sk,K,M A (z)[m]( u ) = SkO + SklU' 1 + Sk2U~ 2 H 

for suitable € End {M.\_{z)[m\). 

Theorem 8.2. The following statements hold. 

(i) The operators Sio,S2o, ■■■ ,<S no and Sn,S2i, ■ ■ ■ ,Sni are scalar operators. Moreover, the 
following relations hold: 

N N 

xN + e s *> xN ~ k = n ( x - k & 

k=l i=l 
N N N 

II (—«■') Ej^:- 

fc=l i=l j=l J 

(mJ .For fe = 1, . . . , JV — 1, w;e have 

n 1 

«5 fe)Q) M A ( z )(«) = 5 fc (u) J] ( ^_^ )fc , (8.2) 

where Sk(u) is a polynomial in u of degree nk. Moreover, the operators 
n 1 

Sk,r = Sk{z r ) -, r-r 

are scalar and 

N N-r-1 N 

(-i) fc s k , s n ^ - j) = n ( d - A i' ] -N+i). (8.3) 

fc=0 J=0 i=l 

Proof. Part (i) follows from Proposition B.l in [6]. The existence of presentation (|8.2[) follows 
from the definition of the universal differential operator. To prove equation (|8.3p it is enough 
to notice that the leading singular term of the universal differential operator at u = z r is equal 
to the leading singular term of the universal differential operator associated with one evaluation 
module M A ( r )(z r ), which in its turn expresses via the quantum determinant. ■ 
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8.3.1. Assume that K is the zero matrix. Then the associated transfer matrices preserve 
Sing A4A( ;Z )[ m ] and we may consider the universal differential operator 

N 

^ K=0, Sing M A (z)[m](. u ,9 u ) = ^ (~l) k S k (u) d^~ k 

k=0 

acting on Sing .Ma valued functions of u. 

Theorem 8.3. For k = 1, . . . , N, the coefficients Sk(u) have the following Laurent expansion 
at u = oo: 

S k {u) = S kfi u~ k + Sk^- 1 + ■■■ , 

where the operators Si^, . . . ,<Sjvo are scalar operators. Moreover, 

N N-k-1 N 

E(-l)\oII (d-j) = H(d-m s -N + s). 

k=0 j=0 s=l 

The proof of Theorem 18,31 is similar to the proof of Theorem 13.71 
8.4 First main result in the Gaudin case 

Theorem 8.4. Let K be an N x N -matrix and M\(z) the tensor product of polynomial 
Ql N [x]-modules. Consider the universal differential operator T^K,MK{z){ u ^d u ) associated with K 
and M\(z). Then the kernel ofT>xj^ A ^(u,d u ) is generated by quasi- exponentials. 

A statement of this type was conjectured in [1]. Theorem 18.41 will be proved in Section [TUl 

9 Continuity principle for differential operators 
with quasi- exponential kernel 

9.1 Quasi-exponentials 

Let V be a complex vector space of dimension d. Let Aq(u), . . . , An(u) be End (V)-vahied 
rational functions in u. Assume that each of these functions has limit as u — ► oo and Aq(u) = 1 
in End (V). Then the differential operator 

k=0 

acting on ^/-valued functions in u, will be called admissible at infinity. 

For every k, let A k {u) = A^q + A^u" 1 + A^u" 2 + • • • be the Laurent expansion at infinity. 
Consider the algebraic equation 

det (x N + x N ' l A lfl + ■■■ + xAtv-1,0 + ^tv,o) = (9.1) 

with respect to variable x. 

Lemma 9.1. If a nonzero V -valued quasi- exponential e Xu (u d Vd + u d ~ 1 v ( i-i + • • • + vq) lies in 
the kernel of an admissible at infinity differential operator V, then X is a root of equation (\9.l\) . 
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9.2 Continuity principle 

Let Ao(u,e), . . . ,Aj^(u,e) be End (F)-valued rational functions in u analytically depending on 
e € [0, 1). Assume that 

N 

• for every e 6 [0, 1) the difference operator T> e = ^ A^{u, e) d^~ k is admissible at infinity, 

fc=o 

• for every e £ (0, 1) the kernel of D e is generated by quasi-exponentials, 

• there exists a natural number m such that for every e € (0, 1) all quasi-exponentials 
generating the kernel of T> 6 are of degree less than m. 

Theorem 9.2. Under these conditions the kernel of the differential operator P e =o is generated 
by quasi-exponentials. 

The proof is similar to the proof of Theorem 14.41 

10 Bethe ansatz in the Gaudin case 

10.1 Preliminaries 

Consider as the gt^-module with highest weight (1, 0, . . . , 0). For complex numbers z\, . . . , 
z n , denote 

M(z) =C N ( Zl )® ■■■®C N (z n ), 

which is the tensor product of polynomial fll/vM -modules. Let 

M(z) = ® mi ^...^m N ^oM(z)[m] 

be its q tTv-weight decomposition with respect to the Cartan subalgebra of diagonal matrices. 
Here m = (mi, . . . , mjv). 

Assume that K = diag (K\, . . . , Kn) is a diagonal N x iV-matrix with distinct coordinates 
and consider the universal differential operator 

N 

V K ,M{ z ){^d u ) = Y J {-^) k S KK ,M{ z )(u)d^- k . 

k=0 

associated with M.(z) and K. Acting on .M(z)-valued functions, the operator T> K M ^{u,d u ) 
preserves the weight decomposition. 

In this section we shall study the kernel of this operator restricted to A4 (z) [m]- valued func- 
tions. 

10.2 Bethe ansatz equations associated with a weight subspace 

Consider a nonzero weight subspace M(z)[m]. Introduce I = . . . , Zjv-i) with lj = mj + i + 
• • • + 7717V- We have n ^ l\ ^ • • ■ ^ In-i ^ 0. Set lo = In = and I = l\ + ■ ■ • + ijv-i- In what 
follows we shall consider functions of I variables 

f _ (Jl) Al) M ,(2) (N-l) (iV-l)x 

L — \ L \ ! l l ! • • ■ 5 %}■■■■> b \ ! • • ■ 5 i l N _ 1 ) ■ 

The following system of I algebraic equations with respect to I variables t is called the Bethe 
ansatz equations associated with jM(^)[m] and matrix K, 

n - h ^ h 2 

«=1 l 3 Zs J'=l l j V J'=l l j V 



Generating Operator of XXX or Gaudin Transfer Matrices Has Quasi-Exponential Kernel 27 



la + 1 la — 1 la ^ 

Yl -joj a^+T) + Z T^iy ~ Z >y = Ka+1 ~ Ka ' 



ijV-2 j JjV-1 2 



Z f (JV-2) N K N-1- 

r i '/ i'=i './ V 

Here the equations of the first group are labeled by j = 1, . . . , l±, the equations of the second 
group are labeled by a = 2, . . . , N — 2, j = 1, . . . , l a , the equations of the third group are labeled 
by j = 1, ... , Zat-i- 

10.3 Weight function and Bethe ansatz theorem 

We denote by uj(t,z) the universal weight function of the Gaudin type associated with the 
weight subspace A4(z)[m]. The universal weight function of the Gaudin type is defined in |12j . 
A convenient formula for uj(t, z) is given in Appendix in [TTj and Theorems 6.3 and 6.5 in [llj . 

If t is a solution of the Bethe ansatz equations, then the vector oj(t, z) € ,M(;s)[m] is called 
the Bethe vector associated with t. 

Theorem 10.1. Let K be a diagonal matrix and t a solution of the Bethe ansatz equations 
(jlO.ip . Assume that the Bethe vector u(t, z) is nonzero. Then the Bethe vector is an eigenvector 
of all transfer-matrices Sk,K,M(z)(, u )> k = 0, . . . , N . 

The statement follows from Theorem 9.2 in [6]. For k = 1, the result is established in |10| . 

The eigenvalues of the Bethe vector are given by the following construction. Set 

xHu,t,z) = K l + ±^-J2-^, 

8=1 i=l u ~ l i 

la— I ^ la 

x a (u, t,z)=K a + J2 T^T) " E 75) > 

i=i u — t i i= i u — t i 

for a = 2, . . . , N. Define the functions Afc(ii, t, z) by the formula 

N 

(d u - X \u, t, z))... (d u - X N (u, t, x)) = E ("l) fe h(u, t, z) d^~ k . 

k=0 

Then for k = 0, . . . , N, 

s k,Q,M{z)( u ) w (*> z ) = h(u,t,x)LL>(t,x), 
see Theorem 9.2 in [6]. 

10.4 Differential operator associated with a solution 
of the Bethe ansatz equations 

Let t be a solution of system (jlO.ip . The scalar difference operator 

N 

V- t ( U ,T) =E(" 1 ) fc h(u,t,z)d^- k 
k=0 

will be called the associated fundamental differential operator. 

Theorem 10.2. The kernel of T>i(u,t) is generated by quasi- exponentials of degree bounded 
from above by a function in n and N . 

This is Proposition 6.4 in [9], which is a generalization of Lemma 5.6 in [TJ. 
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10.5 Completeness of the Bethe ansatz 

Theorem 10.3. Let zi,...,z n and K = diag (Ki, . . . , Kn) be generic. Then the Bethe vectors 
form a basis in A4(z)[m]. 

Theorem 110.31 will be proved in Section 110.61 

Corollary 10.4. Theorems 110.31 and 19.21 imply Theorem[ 



Proof. Theorems 110.21 and 110.31 imply that the statement of Theorem 18.41 holds if the tensor 
product M{z) is considered for generic z and generic diagonal K. Then according to the remark 
in Section 8.2.3, the statement of Theorem 13.81 holds if the tensor product M.(z) is considered 
for generic z and generic (not necessarily diagonal) K. Then the remark in Section 8.1.3 and 
Theorem 19.21 imply that the statement of Theorem 18.41 holds for the tensor product of any 
polynomial finite-dimensional fllyvM -modules and any K. ■ 

10.6 Proof of Theorem flTOI 

10.6.1. For a nonzero weight subspace A^(z)[m] C M.{z) denote by d[m] its dimension. Let 
n ^ h ^ ■ ■ ■ ^ Zat-i ^ be the numbers defined in Section [10. 21 

A vector a = (ai,...,a n ) with coordinates a, from the set {0, 2, 3, . . . , iV} will be called 
admissible if for any j = 1, . . . , N — 1 we have lj = a, | i = 1, . . . , n, and ai > j }. In other 
words, a is admissible if rrij = #{ a,{ \ i = 1, . . . , n, and at = j}. 

If a = (oi, . . . , a n ) is admissible, then for j = 1, . . . , N — 1, there exists a unique increasing 
map p at i : {1, . . . , k} -> {1, . . . , n} such that #p~^(j) = 1 if % > i and #/0~](j) = if aj < i. 

Let v = (1,0, ... ,0) € be the highest weight vector. The set of vectors e a v = e ai ^v (£> 
■ ■ ■ <8> ea n ,iv £ -M(z), labeled by admissible indices a, form a basis of M.(z)[m\. 

10.6.2. Let z\, . . . , z n be distinct numbers. Assume that K = diag (1/q, 2/q, . . . , N/q), where q 
is a small nonzero parameter. Then the Bethe ansatz equations (llO.ip take the form 

n 1 12 1 ll 2 1 

Eli) ~ + (i) (2) " Yl 71) 7iy = o' (10 ' 2) 
,=i - z s /=1 - f=1 t j - t jt « 



El ^ — "\ 1 ^ — "\ L 1 

7S) (a+i) + 2^ IS) + (o-i) ~ 2^ T5o >) = o' 

2^ (7V-1) (7V-2) 2-> (AT-1) (7V-1) Q ' 

Lemma 10.5. For any admissible index a = (ai, . . . , ajv) a^d small nonzero q, there exists 
a solution i(a,q) = (t^\a,q)) of system (| 1 . 2 1) suc/i £/ia£ 

f (a, 5) = + £ - g + 0(, 2 ) (10.3) 

£, = 1 Pa.iU) 

for every i, j. 

Lemma 10.6. As q tends to zero, the Bethe vector u(t (a, q), z) has the following asymptotics: 

u(t(a,q),z) = C a q~ l e a v + 0(q~ l+1 ), 
where C„ is a nonzero number. 
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The lemma follows from formula (|10.3j) and the formula for the universal weight function 
in[TI]. 

Lemma 110.61 implies Theorem 110.31 

11 Comparison theorem in the Gaudin case 

11.1 The automorphism Xf an d the universal differential operator 

oo 

For a series f(u) = f s u~ s ~ l with coefficients in C we have an automorphism 

s=0 

Xf ■ &n[x] ~> fllivW, 4? ^ e { a l } + f s 5 ab . 

For a gljvfxj-module V, denote by V? the representation of fll^M on the same vector space 
given by the rule 

X\y f = (Xf(X))\ V 

for any X 6 fltjvW- 

Lemma 11.1. Let A = (Ai, . . . , Ajy) be a dominant integral gljy -weight and a a complex number. 
Denote A(a) = (Ai + a, . . . , Ajy + a). Let V = M\(z) be an evaluation module and f(u) = 
Then Vf = M A{a) (z). 

11.1.1. Let A = (A^ 1 ), . . . , A( n )) be a collection of integral dominant gl^y-weights, where A® = 
(Aj , ...,A^) for i = l,...,n. Let a = (ai,...,a n ) be a collection on complex numbers. 
Introduce the new collection A(a) = (A( 1 )(a 1 ), . . . , A^(a n )). 

For complex numbers z\, ■ ■ ■ , z n , consider the tensor products of evaluation modules 

Ma(z) = M AW ( Zl ) ® ••• ® M A{n) (z n ), 
M A(a) (z) = M A( i )(ai) («i) ® • • • ® Af A(n)(an) (z n ). 

Theorem 11.2. For any matrix K , an M\(z)-valued function g(u) belongs to the kernel of the 
differential operator L)x : Ma(z)( u j9 u ) if and only if the function C(u)g(u) belongs to the kernel 
of the differential operator T^k^ ^(u,d u ), where 

n 

C(u) = H( U -z i ) ai . 

i=i 

11.1.2. Assume that Ma(z) is the tensor product of polynomial ^[^[xj-modules and a = 
(ai, • • • , ajv) are non-negative integers. Then C(u) is a polynomial and M-A{a){ z ) is the tensor 
product of polynomial gt^r [^-modules. In this case Theorem 111.21 allows us to compare the 
quasi-exponential kernels of the differential operators _A/f A 0)( n i ®u) anci ^k,m a ^(z)( u ^ d u ). 

12 The kernel in the tensor product of evaluation modules 
in the Gaudin case 

12.1 The second main result in the Gaudin case 

Let K = diag (K\, . . . , Kn) be a diagonal matrix. Let A = (A^ 1 ), . . . , A( n )) be a collection of 
integral dominant gl^-weights, where A^ = (A^, . . . , A®) and A^J = for i = 1, . . . , n. 

Choose a weight subspace A / [A(- Z )[ m ] C M\(z) and consider the universal differential oper- 
ator T> K ^ /tA r z \t m ](u,d u ) associated with K and M\(z)[m]. 
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Theorem 12.1. Assume that the numbers K\, . . . , Kn are distinct. Then for any i = 1, . . . , N 

and any nonzero vector vq £ M\(z)[m], there exists an A4\(z)[m]-valued quasi- exponential 

e KiU (v u mi + vm^- 1 + --- + v mi ) (12.1) 

which lies in the kernel of T>x_^M A {z)\m}{ u )d u ). Moreover, the C-span of all such quasi- exponen- 
tials is the kernel of ' ^ir, jvf A ( z )[rn.] (it, <9 U )- 

Proof. On one hand, by Theorem 18.41 the kernel of I'i^.M aC^OM^' d u ) is generated by quasi- 
exponentials. On the other hand, by Lemma 19.11 and part (i) of Theorem 18.21 any quasi- 
exponential lying in the kernel of T> K M A (z)[m]( u i d u ) must be of the form (|12.ip . where Vq is 
a nonzero vector. Since the kernel is of dimension N ■ dim Ai\(z)[m], there exists a quasi- 
exponential of the form (|12.ip with an arbitrary nonzero vq. ■ 

Theorem 12.2. For i = 1, . . . , n, j = 1, . . . , N, and any nonzero vector v £ M\(z)[m], there 
exists an M. \(z)[m]-valued solution f(u) of the differential equation T> K ^M A {z)[m] (u, d u )f(u) = 
such that 

/(«) =v(u- z^ )+N ~> + 0((« - Zi )^ +N ^ +1 ). 

The proof easily follows from part (ii) of Theorem 18.21 

12.1.1. Assume that K is the zero matrix. Consider the subspace Sing A^a^)^] cA4\(z)[m] 
of gljy-singular vectors and the associated universal differential operator TJ K=0 SingMA ( z ^ m }(u,d u ). 

Theorem 12.3. For any i = 1, . . . ,N and any nonzero vector vq £ Sing M\(z)[mi, . . . , m^}, 
there exists a Sing Ai\(z)[m]-valued polynomial 

VQ u m l+ N- t + ^ u m l+ N-i-l + . . . + Wjv _. (12 . 2 ) 

which lies in the kernel ofD K=0SingMA ^ z )[ m j(u,d u ). Moreover, the C-span of all such polyno- 
mials is the kernel of £># = o,Smg M A (z)[m](% 9 U ). 

Proof. On one hand, by Theorem 18.44 the kernel of 'D K=0 ^ n gM A (z)[m]( u , 9 U ) is generated 
by quasi-exponentials. On the other hand, by part (i) of Theorem 18. 2\ any quasi-exponential 
lying in the kernel of T> K=0 Sing M A (z)[m] (u, d u ) must a polynomial. By Theorem 18.31 such a 
polynomial has to be of the form indicated in f|12.2|) . where vq is a nonzero vector. Since the 
kernel is of dimension iV-dim Sing A4\(z)[m], there exists a polynomial of the form (|12.2p with 
an arbitrary nonzero vq. I 

12.2 Kernel of the fundamental differential operator associated 
with an eigenvector 

Assume that v £ Ai\(z)[m] is an eigenvector of all transfer matrices, 

Sk,Ma(z)[™] ( u ) v = h,v(, u ) v > k = 0, . . . ,N. 
Then the scalar differential operator 

TV 

iN-k 



v v { u ,d u ) = Y J {-if\MdX 



k=0 



will be called the fundamental differential operator associated with the eigenvector v. 

Theorems 112. 11 [12. 21 and 112.31 give us information on the kernel of the fundamental differential 
operator. 
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Corollary 12.4. 

(i) Assume that K is diagonal with distinct diagonal entries. Then the kernel of T> v (u,d u ) is 
generated by quasi- exponentials e KlU pi(u), . . . ,e KNU p^{u), where for every i the polyno- 
mial Pi{u) £ C[u] is of degree mi. 

(ii) Assume that K is the zero matrix and the eigenvector v belongs to Sing A4\(z)[m]. Then 
the kernel ofT> v (u, d u ) is generated by suitable polynomials pi(u), . . . ,pn( u ) °f degree mi + 
N — 1, . . . ,m,N, respectively. 

Corollary 12.5. For i = 1, . . . ,n, j = 1, . . . , N , there exists a solution f(u) of the differential 
equation T> v (u,d u ) f(u) = such that 

/(«) = (ti - Zi)*?*"-' + 0((u- *i) A " + ^ +1 ). 
If v is a Bethe eigenvector, then these two corollaries are proved in [7] and [9]. 
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